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Abstract. This is an overview about natural sample spaces for differential equations driven 
by various noises. Appropriate sample spaces arc needed in order to facilitate a random dy- 
namical systems approach for stochastic differential equations. The noise could be white 
or colored, Gaussian or non-Gaussian, Markov or non-Markov, and semimartingale or non- 
semimartingale. Typical noises are defined in terms of Brownian motion, Levy motion and 
fractional Brownian motion. In each of these cases, a canonical sample space with an appro- 
priate metric (or topology that gives convergence concept) is introduced. Basic properties 
of canonical sample spaces, such as separability and completeness, are then discussed. 

Moreover, a flow defined by shifts, is introduced on these canonical sample spaces. This 
flow has an invariant measure which is the probability distribution for Brownian motion, or 
Levy motion or fractional Brownian motion. Thus canonical sample spaces are much richer 
in mathematical structures than the usual sample spaces in probability theory, as they have 
metric or topological structures, together with a shift flow (or driving flow) defined on it. 
This facilitates dynamical systems approaches for studying stochastic differential equations. 



1. Random dynamical systems 

Stochastic differential equations (SDEs) or stochastic partial differential equations (SPDEs) 
arise as mathematical models for complex systems under various random influences in en- 
gineering and science. Here we only consider random dynamical systems defined by SDEs. 
Such SDEs define random dynamical systems (RDS) with appropriate sample spaces, much 
as ordinary differential equations define deterministic dynamical systems. 
Theory of random dynamical systems allows to discuss the qualitative behavior of stochastic 
systems that are not only driven by a white noise, Markov processes and semimartingales, but 
also driven by non-Markov processes or by non- semimartingales (e.g., fractional Brownian 
motion). To analyze these more general noise cases, appropriate sample spaces and ergodic 
theory play an important role. In this article, we discuss canonical or natural sample spaces 
for SDEs with various noises. 

We recall the definition of a random dynamical system (RDS) in the state space H = R n , 
with the underlying probability space (Q, J 7 , P), and with time t varying in T = R = 
(—00, 00) or T = R + = [0, 00), as in Arnold pQ. The state space R n is equipped with the 
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Euclidean norm (or length) \x\ = \J x\ + ■ • • + x 2 n and the usual scalar product < x,y >= 
X\V\ H V x n y n . 

Note that a deterministic dynamical system on the state space if is a mapping if) : T x H —> 
H, (t, x) i — > ifj(t, x), such that the flow property is satisfied: 

ip(0,x)—x, ip(t + s,x) = ip(t,ip(s,x)), 

for all t, s G T and x G H. 

For a random dynamical system, we need an extra ingredient, namely, a model for the noise. 
Moreover, the flow property has a twist (thus called cocycle property) due to the effect of 
noise. 

Definition 1.1. (Random dynamical system) 

A random dynamical system (RDS), denoted by if, consists of two ingredients: 

(i) Model for the noise: A driving flow on a probability space (f2,jF, P), i.e., a flow 
(8t)teT on the sample space Q, such that P is invariant, namely 6tP = P for all t G T, and 
(t, uj) i— > 9 t uj is measurable from TxO to Q. 

(ii) Model for the evolution: A cocycle if over 9, i.e. a measurable mapping (p : T x Q x 
H — > H, (t, uj, x) t— > tp(t, uj, x), such that the family (p(t, uj, ■) = <p(t, uj) : H — > H of random 
mappings satisfies the cocycle property: 

(1) ip(0, uj) = idn, <p(t + s, uj) = ip(t, 9 s uj) o <^(s, uj) for all t, s G T, uj G f2. 

Here the driving flow 9t describes stationary dynamics of noise in an appropriately chosen 
sample space (see below). The mathematical model for noises in engineering and science is 
usually a stationary generalized stochastic process [22] . 

When (t, x) i— > tp(t,uj,x) is continuous for all uj G f2, we say that tp is a continuous RDS. 
Since the continuity in space x is quite common for RDS generated by stochastic differential 
equations, we usually do not specifically mention this spatial continuity. We often call <p a 
continuous-time or discrete-time RDS when it is continuous or discrete in time t. 
It follows from [l] that ip(t,uj), t G R, is a homeomorphism of H and 

(p(t,uy l = <f(-t,9 t uj). 

2. Dynamical systems driven by white noises 

2.1. Brownian Motion. The physical phenomenon Brownian motion^ is due to the inces- 
sant hitting of pollen by the much smaller molecules of the liquid. The hits occur a large 
number of times in any small time interval, independently of each other and the effect of a 
particular hit is small compared to the total effect [13] . The physical theory of this motion, 
set up by Albert Einstein in 1905, suggests that the motion is random, and has the following 
properties: 

i) the motion is continuous; 

ii) it has independent increments; 



The phenomenon was first observed by Jan Ingenhouz in 1785, but was subsequently rediscovered by 
Brown in 1828, according to sources used by Eric Weisstein's World of Physics, which can be found on the 
Internet at |http:/ /scienceworld.wolfram.com/physics/BrownianMotion.html 
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A Scalar Brownian Motion 
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FIGURE 1. A sample path of Brownian motion B(t) 



iii) the increments are stationary and Gaussian random variables. 
Figure 1 shows a sample path of the Brownian motion. 

Intuitively speaking, property i) says that the sample path of the Brownian motion is con- 
tinuous. Property ii) means that the displacements of a pollen particle over disjoint time 
intervals are independent random variables. Property iii) is natural considering the Central 
Limit Theorem. 

We now describe the Brownian motion in the mathematical language, i.e., introduce the 
first definition of the Brownian motion[22j [2j [6]. 

Definition A: A stochastic process {B t (u) : t > 0} defined on a probability space J 7 , P) 
is called a Brownian motion or a Wiener process if the following conditions hold: 

1) B (u) = a.s.; 

2) the sample paths t — > B t {oj) are a.s. continuous; 

3) B t {uj) has stationary independent increments; 

4) the increments B t (u) — B s (uj) has the normal distribution with mean and variance t — s, 
i.e. B t (uo) - B s (uo) ~ N(0, t - s) for any < s < t. 

Since the stochastic process B t (uj) is a mapping from a probability space to a metric space 
and is governed by its law, i.e. the probability measure on the metric space induced by B t (u), 
we want to find it in order to have a better understanding of the Brownian motion. In fact, 
we can find the finite dimensional distribution of the Brownian motion using condition 3) 
and 4) in Definition A, and this gives another definition of the Brownian motion as we shall 
see. We first write down the second definition [6] and then prove it is equivalent to Definition 
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A. 

Definition B: A stochastic process {B t (u) : t > 0} defined on a probability space (Q, J 7 , P) 
is called a Brownian motion or a Wiener process if the following conditions hold: 
1') B (u) = a.s.; 

2') the sample paths t — > B t {oj) are a.s. continuous; 

3') for any finite sequence of times = to < t\ < ti < ■ ■ ■ < t n and Borel sets B±, ■ ■ ■ ,B„cR 



P{B tl (cu)EB 1 ,--- ,B tn (u)EB n } 
p(ti, 0, xx)p(t 2 - tx,xt, x 2 ) ■ ■ ■ 



(2) • • -p(t n - t n _ x ,x n -x,x n )dxi ■■■dx r 



where 



(3) p(t;x,y) 



defined for any x, y G M and t > is called the transition density. 

We can see the conditions l')>2') in Definition B are completely the same as 1),2) in Defini- 
tion A, respectively, the only thing we need to do in the proof of the equivalence of the two 
is to show 3), 4) & 3'). 

Proof 3),4)=> 3') 

Firstly, we show that Brownian motion B t has Markov property, by proving that the con- 
ditional distribution of B t + S given JF 4 is the same as that given B t , in terms of moment 
generating function [10J. In fact, 



E(e uBt+ °\F t ) = E(e u[{Bt+s - Bt)+Bt] \F t ) 
= e uBt E(e u{Bt+s ~ Bt \Tt) 
= e uBt E(e u( - Bt+a ~ Bt ^) 

_ e uBt e u 2 s/2 

= e uBt E{e u{Bt+a - Bt) \B t ) 
= E(e uBt+s \Bt). 
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Secondly, we compute the joint distribution of the Brownian motion. 

P{B{t k+l ) <x k+1 ,B(t k ) < x k } 
= P{B(t k ) < x k , [B{t k+1 ) - B(t k )} + B(t k ) < x k+1 } 
= P{B{t k ) < x k , [B(t k+1 ) - B(t k )} < x k+1 - B{t k )} 

= P{B{t k+l ) - B{t k ) < x k+l - B(t k )\B(t k ) < x k }P{B(t k ) < x k } (Conditional probability) 

X k t-X k+1 -X ^ y 2 ^ ^2 

— p ^k+i-tk) dy p~2*kd/r 

^2ir(t k+1 - t k ) V^h 



oo J — oo 



x k ^ x 2 fXk+l ^ (y-x) 2 

e~^dx / e ^k+i-'k) dy 

oo V2vrt fc ,/„ 00 y/2n(t k+1 - t k ) 

%k r^k+i 

p(t k ;0,x)dx p(t k+1 -t k ;x,y)dy. 

oo J — oo 

Finally, we have the finite dimensional distribution. 

P{B(h) < x u B(t 2 ) < x 2 , ■ ■ ■ , B{t n ) < x n } 
= P{B{h) < x x \B{t Q ) = 0}P{B(t 2 ) < x 2 \B(t 1 ) < x x } 

P{B{h) < XslBih) < x 1; B{t 2 ) < x 2 } ■ ■ ■ P{B{t n ) < x^Bih) < x x - - ■ fl(t n _i) < x n ^} 
= P{B{h) < x x \B{t Q ) = 0}P{B(t 2 ) < x 2 \B{U) < x x } 

P{B{h) < x 3 \B(t 2 ) < x 2 } ■ ■■P{B(t n ) < x n \B{t n _y) < x n ^} (Markov property) 

PjBjh) < xjPjBjh) < x 2 , gfti) < xx} ■ ■ ■ P{B(t n ) < x n , ggg-i) < x n ^} 
P{B(h) < Xl }P{B(t 2 ) < x 2 } ■ ■ ■ P{B(t n ^) < x„_x} 

pX2 fin 

p(ti;0,yi)dyi p(h -t 1 ;y ll y 2 )dy 2 - ■ ■ p(t n - t n _i;y n _ u y n )dy n 



which is obviously the same as 3') . □ 
Conversely, we need to show 3')=^ 3),4) 

However, this part of work is completely done in [6] (see page 153-155). 



2.2. Wiener Measure. In Section 1, we treat the Brownian motion as a stochastic pro- 
cess, i.e., a collection of time-parameterized random variables. Since a stochastic process £ 
is governed by its law, i.e., the probability measure /x := P^~ l on the space it maps to, one 
may ask the question why the Brownian motion can be determined only by its finite dimen- 
sional distributions although we have proved the so-defined Bownian motion (Definition B) 
coincides with the definition describing the phenomenon (Definition A). Motivated by this 
question, we shall treat the Brownian motion as a "random variable" , which we call random 
function [TT], and this point of view introduces the Wiener measure, the probability mea- 
sure defined on an appropriate space which gives us the Brownian motion. We will sketch 
the ideas showing the existence and uniqueness of this special probability measure, starting 
from three different spaces, i.e, R[0, oo), C[0, oo), C*o[0, oo), the spaces of arbitrary functions, 
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continuous functions, and continuous functions passing zero at time 0, defined on [0, oo), 
respectively. 

First approach This way of showing the existence and uniqueness of the Wiener mea- 
sure is the one most often used to construct a Markov Process and is rather technical. The 
main idea is the following: First define a set function on the algebra generated by the cylinder 
sets in M[0, oo) according to the finite-dimensional distribution of Brownian motion. Note 
that it is a set function rather than a probability measure since it is just finitely additive 
(not countably additive); then by the celebrated Kolmogorov extension theorem [22], this set 
function is uniquely extended to a probability measure on the a— algebra generated by the 
algebra mentioned above. 

More precisely, we now give the definitions and theorems. 

• K[0, oo) denotes the set of all real- valued functions on [0, oo). 

• Cylinder set (22] A subset of R[0, oo) of the form 

(4) A = {uo e K[0, oo) : (ufa), . . . , u{t n )) e B n } 

where B n is a Borel subset of WL n is called a cylinder set. 

Fixing tx,...t n but varying B n over the entire Borel subsets of K n , the class of such all 
cylinder sets forms a a— algebra Q3** 1 ' "'*™*'. 

• Set function A set function $t 1 ,...,t n is defined on the measurable space (R[0, oo), QS^ 1 ''"'*™)), 
given by 

$ tl ,... A (A) :=P{(B h ,...,B tn )eB n } 

By varying the choice of the finite time points {t\, . . . , t n } C [0, oo), we get a class of such 
set functions <3> := {^ti,. ..,*„}, and this class is independent of the cylindrical expression of 
the functions in R[0, oo), i.e. if the cylinder set A of (Jlj) has another expression, say 

A = {u e K[0, oo) : (w(s0, . . . , w(s m )) G B m }, 

then 

$t lt ...,t n {A) = 'I'... ,,. 

• Algebra and er-algebra 

■ it[0, oo) denotes the algebra of subsets of R[0, oo) consisting all cylinder sets 

■ Q5[0, oo) denotes the smallest a— algebra containing lt[0, oo) 

■ $ is a finitely additive measure on (M[0, oo),il[0, oo)) such that the restriction of $ to 
«B(*i'-'*») coincides with $ tl t 

Theorem 2.1 (Kolmogorov extension theorem). The set function $ on (M[0, oo),it[0, oo)) 
is uniquely extendable to a probability measure $ on (R[0, oo), 23[0, oo)). 
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It has been proved [22], [32] that there exist a unique probability measure P on (R[0, oo), 03 [0, oo)), 
under which the coordinate mapping process 



satisfies condition 3) and 4) of Definition A in section 2.1, so it does not introduce a "Brow- 
nian motion" as we defined. Fortunately, by another famous theorem of Kolmogorov|41j . 
the continuity problem has been solved. 

Theorem 2.2 (Kolmogorov continuity theorem). Suppose that the process X = X{t} t >o 
satisfies the following condition: for allT > there exist positive constants a, (3,D such that 



then there exists a continuous modification of X. 

Modification |4T] Suppose that {Xt}, {Y t } are stochastic processes on (Jl,!F,P), then we 
say that {X t } is a modification of {Y t } if 



Note that if X t is a modification of Y t} then they have the same finite-dimensional distribu- 
tions. 

Now there is one problem needs to be solved: Why is B (u) = u(0) = a.s.? 

Definition [11] Wiener measure, denoted by fiw, is a probability measure on a measurable 
space (C, C) having the following two properties: 

i) each ui{t) £ C is normally distributed under fi w with mean and variance t, i.e, 



For t = this is interpreted to mean that fi w {u(0) = 0} = 1. 

ii) the stochastic process {uj(t) : t > 0} has independent increments under fi w , i.e., V0 < 
to <t\ < ■ ■ ■ <t n , Lo(ti) —u)(to),u(t2) —co(ti), ■ ■ -uj(t n ) —uj(t n -i) are independent under /i w . 

Remark It may not seem so obvious that this approach in fact proves the existence of the 
Wiener measure. Also, we cannot obtain the Wiener measure simply by assigning measure 
one to C[0, oo); see [32]. One might hope to construct the measure directly on C[0, oo). 
Indeed, this is the main idea of the second approach we are going to present. 

Second approach There are some advantages if we start from C[0, oo) since we can make 
it Polish (complete and separable) by assigning an appropriate metric. It has been proved 
that the Borel a— algebra generated by the open sets in C[0, oo) is equal to the a— algebra 
generated by all the cylinder sets [32]. Thus there is a totally different way of constructing 
Wiener measure. 



B t (u) := u(t); u £ R[0, oo), t > 



E[\X t -X B \ a ] <D - |t-s| 1+/3 ; 0< s,t<T 



P({u;X t (u)=Y t (u)}) = lVt 
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Metric on C[0,oo) 



(5) 



p(u a ,uj b ) := ^ 



2 0<t<n 



n=l 



This metric introduces the topology of uniform convergence on compact intervals. For con- 
venience, we denote C[0, oo) by C and its Borel cr-algebra by C. 

The main idea is to construct a sequence of probability measure {P n } on (C,C) such that 
Pn =^ f^w Since ([2]) gives the f.d.d of the Wiener measure fi w , we must let the f.d.d of {P n } 
weakly converges to those of /j, w . Although weak convergence in C need not follow the weak 
convergence of the f.d.d alone in general, it does if we add the condition that {P n } is tight 
[3"2T [TTj . In fact, if a metric space is separable and complete, then each probability measure 
on the measurable space is tight (see Theorem 1.3 in [TT]). It can be shown that equipped 
with the metric defined by (j5J), the canonical space is separable and complete [TT]. For detail, 
see [32] [UJ. 

Remark What is still not natural is that Bq(uj) = u)(0) = a.s. So next we will talk 
about constructing Wiener measure on Cq[0, oo). It appears that the above two different 
approaches may be adapted to this case. 

We can define the Brownian motion B t for t G R as follows: Taking two independent 
Brownian motions B t and B t , we define 



We can work on the space Co(R, R n ) for two-sided Brownian motion. The Wiener measure 
defined above should be similar defined in this space pQ. 

2.3. Canonical sample space. We consider a SDE 



The canonical sample space is Q := Co(R'R n ), space of continuous functions that are zero at 
time zero, equipped with the compact open topology, the Borel a— field T := £>(Co(R + , R n )), 
and Wiener (probability) measure fiw 

We introduce a driving flow 9 t on this canonical sample space Q is given by the Wiener shifts 

9 t oj{-) = oj(- + t) -u(t), teR, u 6 n = C Q (R,R d ). 
In this case the measure fiw is invariant [TJ, i.e., 

fi w {e;\A)) = fi W {A) 

for all A G T. Moreover, this invariant measure is actually ergodic with respect to the flow 



(6) 




(7) 



dX t = b(X t )dt + a(X t )dB(t). 



em- 
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3. Dynamical systems driven by colored noises 

Colored noise, or noise with non-zero correlation ('memory') in time, are common in the 
physical, biological and engineering sciences [2D]. A good candidate for modeling colored 
noise is the fractional Brownian motion. 

3.1. Fractional Brownian motion. A fractional Brownian motion (fBM) B H (t), i 6 R, 
with H G (0, 1) the Hurst parameter, is still a Gaussian process. But it is characterized 
by the stationarity of its increments and a memory property. The increments of the frac- 
tional Brownian motion are not independent, except in the standard Brownian motion case 
(H = |). Thus it is not a Markov process except when H — \. Specifically, B H (0) = 
a.s., mean EB H (t) = 0, covariance E[B H (t)B H \s)\ = l(\t\ 2H + \s\ 2H - \t - s\ 2H ), and vari- 
ance Var [B H (t) - B H (s)] = \t - s\ 2H . It also exhibits power scaling and path regularity 
properties with Holder parameter H, which are very distinct from Brownian motion. The 
standard Brownian motion is a special fBM with H = 1/2. Figure 2 is a sample path of the 
fractional Brownian motion with H = 0.25. 




t 



FIGURE 2. A sample path of fractional Brownian motion B (t), with H = 0.25 

Note that fBM B H (t) is non-Markov and non-semimartingale. Thus the usual stochastic 
integration [45] is not applicable, and other integration concepts are needed [40, [39] . 

3.2. Canonical sample space. The stochastic calculus involving fBM is currently being 
developed; see e.g. [iQl [51] and references therein. This will lead to more advances in the 
study of SDEs driven by colored fBM noise: 

(8) dX t = b(X t )dt + a(X t )dB H (t). 

Since the fBM B H (t) is not Markov, the solution process X t is not Markov either. Thus the 
usual techniques from Markov processes will not be applicable to the study of SDEs driven 
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by fBms. However, the random dynamical systems approach, as described in ^T] above, looks 
promising [2S1 EES]- The theory of RDS, developed by Arnold and coworkers [I], describes 
the qualitative behavior of systems of stochastic differential equations in terms of stability, 
Lyapunov exponents, invariant manifolds, and attractors. 

See Theorem 2.3 in [381 [9], the canonical sample space is Q := Co(R, R n ), the set of continu- 
ous functions that is zero at zero, but the probability measure jifBM is generated by Bf under 
the compact-open topology as defined in Section[2l The Borel a— field is T := £>(C*o(R, R n ))- 
We can introduce a flow 9 t on this canonical sample space Q defined by the shifts 

6 t u}(')=u(- + t)-u)(t), teR, uj e n = C (M,M d ). 

In this case the measure [I/bm is invariant, i.e. 

lL fBM {pj\A)) = IM fBM (A) 

for all A e T. 

4. Dynamical systems driven by non-Gaussian noises 

In the last two sections, we considered dynamical systems driven by Gaussian noises (white 
or colored), in terms of Brownian motion or fractional Brownian motion. In this section, we 
discuss differential equation driven by non-Gaussian Levy noises. 

4.1. Levy Motions. Gaussian processes like Brownian motion have been widely used to 
model fluctuations in engineering and science. For a particle in Brownian motion, its sample 
paths are continuous in time almost surely (i.e., no jumps), its mean square displacement 
increases linearly in time (i.e., normal diffusion), and its probability density function decays 
exponentially in space (i.e., light tail or exponential relaxation) [41]. But some complex 
phenomena involve non-Gaussian fluctuations, with peculiar properties such as anomalous 
diffusion (mean square displacement is a nonlinear power law of time) [12] and heavy tail 
(non-exponential relaxation) [57]. For instance, it has been argued that diffusion in a case 
of geophysical turbulence [52] is anomalous. Loosely speaking, the diffusion process consists 
of a series of "pauses", when the particle is trapped by a coherent structure, and "flights" 
or "jumps" or other extreme events, when the particle moves in a jet flow. Moreover, 
anomalous electrical transport properties have been observed in some amorphous materials 
such as insulators, semiconductors and polymers, where transient current is asymptotically 
a power law function of time [4HJ r2T] . Finally, some paleoclimatic data [T5] indicates heavy 
tail distributions and some DNA data [52] shows long range power law decay for spatial 
correlation. 

Levy motions are thought to be appropriate models for non-Gaussian processes with jumps 
[48] . Let us recall that a Levy motion L(t), or L t , is a non-Gaussian process with independent 
and stationary increments, i.e., increments AL(t, At) = L(t + At) — L(t) are stationary 
(therefore AL has no statistical dependence on t) and independent for any non overlapping 
time lags At. Moreover, its sample paths are only continuous in probability, namely, P(\L(t) — 
L{to)\ > 5) as t — > to f° r an y positive 5. With a suitable modification [4], these paths 
may be taken as cadlag, i.e., paths are continuous on the right and have limits on the left. 
This continuity is weaker than the usual continuity in time. 
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This generalizes the Brownian motion B(t) or B t , as B(t) satisfies all these three conditions. 
But Additionally, (i) Almost every sample path of the Brownian motion is continuous in 
time in the usual sense and (ii) Brownian motion's increments are Gaussian distributed. 

Dynamical systems driven by non-Gaussian Levy noises have attracted much attention re- 
cently [U [30, ED]. Under certain conditions, the SDEs driven by Levy motion generate 
stochastic flows [351 E]> an d also generate random dynamical systems (or cocycles) in the 
sense of Arnold pQ. Recently, exit time estimates have been investigated by Imkeller and 
Pavlyukevich [271 I2H] , and Yang and Duan [56] for SDEs driven by Levy motion. This 
shows some qualitatively different dynamical behaviors between SDEs driven by Gaussian 
and non-Gaussian noises. 

Levy motions are named in honor of the French probabilist Paul Levy, who first studied 
them in 1930s. From a mathematical point of view, there are so many reasons why they are 
so important [I], such as: 

• There are many important examples, such as Brownian motion, the Poisson process, stable 
processes, and subordinators. 

• They are generalizations of random walks to continuous time. 

• They are the simplest class of processes whose paths consist of continuous motion inter- 
spersed with jump discontinuities of random size appearing at random times. 

Definition [I] A stochastic process L = (L(t),t > 0) defined on a probability space (Q, J 7 , P) 
is a Levy motion if: 
(LI) L(0) = a.s.; 

(L2) L has independent and stationary increments; 

(L3) L is stochastically continuous, i.e. for all e > and for all s 

\imP(\X(t) - X(s)\ > e) = 

t—*s 

With a suitable modification [3], these paths may be taken as cadlag, i.e., paths are con- 
tinuous on the right and have limits on the left. This continuity is weaker than the usual 
continuity in time. 

Figure 3 is a sample path for a Levy motion. 

One way to understand the structure of the Levy motions is to employ Fourier analysis. It 
may be shown that each L(t) is infinitely divisible. The infinitely divisible random variables 
are characterized completely through their characteristic functions by a beautiful formula, 
established by Paul Levy and A. Ya. Khintchine in the 1930s. The related definitions and 
theorems follows. 

Definition [4] The characteristic function of a stochastic process X(t) taking values in 
M. d is the mapping $ 4 : M. d — > C defined by 

$ t (u) = E(e m - X(t) ) := I e lu - y p t {dy) 

where p t is the distribution of X(t). 
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Figure 3. A sample path for a Levy motion 

Definition [3] X(t) is infinitely divisible if for each n G N, there exists a probability mea- 
sure p t ,n on M. d with characteristic function $ t ,n such that $t(u) = ($t 5 „(M)) n , for each u G M d . 



Theorem 4.1 (The Levy-Khintchine Formula [1]). If L = (L(t),t > 0) is a Levy motion, 
then 



where 
(9) 



1 



3 t?7(u) 



for each t > 0, u G 



?7(m) = ib ■ u — -u ■ au + 



tu-y 



K d -{0} 



1 - I'M"!/ I|y|<l(Z/)]^(^2/) 



/or some 6 G a non-negative definite symmetric d x d matrix a and a Borel measure v, 
called Levy jump measure, on ~R d — {0} for which J Rd _^ y(\y\ 2 A l)z/(rfy) < oo. Here Is is the 
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indicator function of the set S. 

Conversely, given a mapping of the form (Q, we can always construct a Levy motion for 
which = e tr} ^ u ' and call it a Levy motion with the characteristics or generating triple 

(a,b,u). 

Remark: The so-defined Borel measure v in (jH]) is called the Levy jump measure, which 
should not be confused with the probability measure induced by the Levy motion, to be 
defined below. 



The different terms which appear in the Levy-Khintchine formula have a probabilistic signif- 
icance emphasized in [16]. Every Levy motion is obtained as a sum of independent processes 
with three types of characteristics (0, b, 0), (a, 0, 0) and (0, 0, v). Thus, the Levy measure ac- 
counts for the jumps of L and the knowledge of v permits to give a probabilistic construction 
of L; see 



4.2. Canonical sample space. Consider a SDE driven by non-Gaussian Levy noise 
(10) dX t = b(X t )dt + a(X t )dL(t). 

The canonical space has to be enlarged to include all the cadlag functions, i.e. functions 
that are right-continuous and have left limits, defined on R and taking values in M. d . This 
space is denoted as D(FL, R n ). 

We adopt the same point of view as in Section [2] that a stochastic process is also a random 
variable, i.e. 

L t (u) : oo — > D(R, R n ) w -> L(t),t eR. 

Remark 4.2. Here the compact-open metric defined in (jSJ) cannot make -D(R + ,R n ) sepa- 
rable [45]. For example, if f a (t) = l[ a ,oo)(t),fp(t) = 1 [/?,«>)(£)> then p(f a , fp) = 1/2 for all 
a,j3 with < a < (3 < 1, and since there are uncountably many such a,j3, the space is not 
separable. 

However, it can be made complete and separable when endowed with the Skorohod metric 
jl lj . With this special metric, we call D(R, M. d ) a Skorohod space. The Skorohod metric on 
D(R,R d ) is defined as 

oo 1 

d(x,y) := ^(^Cf^.l/I) for all x,y e D 

771=1 

where x m (t) := g m {t)x(t), y m (t) := g m (t)y(t) with 

1, if \t\ < m - 1 



9m(t) ■-- 



m — t, if m — 1 < \t\ < m, 
0, if \t\ > m 



and 



d° m (x,y) := mf 



sup 

-771<S<t<77l 



log 



X(t) - \(s) 



V sup \x(t) - y(X(t))\ 

—m<t<m 
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where A denotes the set of strictly increasing, continuous functions from R to itself. 

The Borel a— field under this topology is denoted as S. For studying the weak convergence 
and tightness in D, the same approach adopted in C can be applied except that the fact the 
natural projections are not continuous need to be noticed [IT] . 

Definition The probability measure, /i^, in (D(R, R n ), S) that makes every element in 
-D(R, R n ) a sample Levy path is called the Levy probability measure. Note that this measure 
is not to be confused with the Levy jump measure v mentioned above. 
We can also introduce a flow 9 = (6t,t G R) on this canonical sample space Q by the shifts 

(11) (6 t uj)(s):=uj(t + s)-uj(t). 

The (Levy) probability measure \xl is invariant under this flow, i.e. 

i u L (6; 1 (A)) = i u L (A) 

for all A G T\ see jl] (Page 325) and [37]. This flow is an ergodic dynamical system [T] with 
respect to the above probability measure /i^. 
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